BASIC ZETA FUNCTIONS AND SOME APPLICATIONS IN 

PHYSICS 

KLAUS KIRSTEN 

o 

£>-». 1. Introduction 

a 

It is the aim of these lectures to introduce some basic zeta functions and their 
uses in the areas of the Casimir effect and Bose-Einstein condensation. A brief 
introduction into these areas is given in the respective sections; for recent mono- 
graphs on these topics see [22j EH [MJ ED ID ES ED E2]- We will consider 
exclusively spectral zeta functions, that is zeta functions arising from the eigen- 
value spectrum of suitable differential operators. Applications like those in number 
theory [3J 2J [531 EH] will n °t be considered in this contribution. 
Q-ij There is a set of technical tools that are at the very heart of understanding 

analytical properties of essentially every spectral zeta function. Those tools are 
introduced in Section 2 using the well-studied examples of the Hurwitz |54| , Epstein 
[38[ [39] and Barnes zeta function [3 [6]. In Section 3 it is explained how these 
different examples can all be thought of as being generated by the same mechanism, 
namely they all result from eigenvalues of suitable (partial) differential operators. 
qq | It is this relation with partial differential operators that provides the motivation 

CO ■ for analyzing the zeta functions considered in these lectures. Motivations come 

CNI | for example from the questions " Can one hear the shape of a drum?" and " What 

l/~) ■ does the Casimir effect know about a boundary?". Finally "What does a Bose 

gas know about its container?" The first two questions are considered in detail in 
Section 4. The last question is examined in Section 5 where we will see how zeta 
functions can be used to analyze the phenomenon of Bosc-Einstcin condensation. 
The Conclusions will point towards recent developments for the analysis of spectral 
zeta functions and their applications. 
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2. Some basic zeta functions 

In this section we will construct analytical continuations of basic zeta functions. 
From these we will determine the meromorphic structure, residues at singular points 
and special function values. 

2.1. Hurwitz zeta function. We start by considering a generalization of the Ric- 
mann zeta function 



(2-i) cum = E ^ 



n 

n=l 



Definition 2.1. Let s <G C and < a < 1. Then for Sfts > 1 the Hurwitz zeta 
function is defined by 



oo .. 



—i(n + a) s 
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Clearly, we have that (h{s, 1) = (r(s). Results for a = 1 + b > 1 follow by 
observing 

oo 
n— 

In order to determine properties of the Hurwitz zeta function, one strategy is to 
express it in term of 'known' zeta functions like the Riemann zeta function. 

Theorem 2.2. For < a < 1 we have 

fc=0 v ' 

Proof. Note that for \z\ < 1 we have the binomial expansion 



(i - z)- = e U^Sz k 



, rww 



So for -fts > 1 we compute 

1 °° 1 i 

( H (s,a) = — + V— . , 

a s *— ' n s (1 + — 



n=l 

oo _. oo 



- 4 + E i 7E(- 1 ) fc Tf4#(- 

n=l fc=o v ' 

= - + V(-l) fc r(s + fc) a fc V^- 
a» ^ J r(s)fc! ^ n s+fc 

fc=0 V ' n=l 

fc=0 v ' 

which is the assertion. □ 

From here it is seen that s = 1 is the only pole of (h(s, a) with Res £#(1, a) = 1. 
In determining certain function values of Ch(s, a) the following polynomials will 
turn out to be useful. 

Definition 2.3. For x <G C we define the Bernoulli polynomials B n (x) by the 
equation 

(2.2) -^—- = Y ^4^z™, where \z\ < 2tt. 
y ' e z - 1 ^ n! ' ' 

Examples are Bq{x) = 1 and -Bi(a;) = a; — 1/2. The numbers B„(0) are called 
Bernoulli numbers and are denoted by B n . Thus 

oo „ 

(2.3) —=— = V -±z n , where |z| < 2tt. 

n=0 

Lemma 2.4. T/ie Bernoulli polynomials satisfy 
(1) 

5 n (x) = E(fc) B ^"" fe ' 
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(2) 

B n {x + 1) - B n (x) = nx"- 1 ifn>l, 

(3) 

(-l) n B n (-x) = B n (x) + nx n '\ 

(4) 

B n (l-x) = (-l) n B n (x). 

Exercise 1. Use relations (|2.2[) and (|2.3|) to show the assertions of Lemma [23 

We now establish elementary properties of ^n{s,a). 

Theorem 2.5. For 5Rs > 1 we have 

oo 

(2.4) c^- y^i__*. 

o 
Furthermore, for k £ Nq we have 

z- , i \ S fc+1 (a) 

Proof. We use the definition of the Gamma-function and have 

oo oo 

,s-l„-u J„, _ \s / +s-l -tA, 



(2.5) r(s) = / u a_1 e _u du = A s / t s ~ x er tx dt. 



This shows the first part of the Theorem, 

oo oo 

OO -. p -. -. oo 

feM = Eri)i rVt(n+a)dt = ri)/ rl S e " i(n+o)rft 

71 U ri ri 71 U 
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rd*. 



r( s ) y i - e - 



Furthermore we have 





1 oo 

1 /\ s - 2 H)e- t ° 1 r._ a H)e-* 



d£+=^- / £ s ~ - — f dt. 



T(s) J e"* - 1 T(s) 7 e- 

o 1 

The integral in the second term is an entire function of s. Given the Gamma- 
function has singularities at s = —k, k G No, only the first term can possibly 
contribute to the properties Cff(~ k, a) considered. We continue and write 

Vis) I £- — ' m\ 



T(s)J e-t-1 T(s)J f- n 



n 



1 ^B n (a) (-1)' 



E 



T(s) ^— ' n! s + ?i — 1 

s 7 n— 
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which provides the analytical continuation of the integral to the complex plane. 
From here we observe again 

ResCff(l,a) = B (a) = l, 

and the second part of the Theorem 

1 B k+1 (a) (-l) fc +! 



( H (-k,a) = lim 



oT(-k + e) (fc + 1)! e 

^o v ' (fc + 1)! e 



Bfc+i(o) 



follows. 



□ 



The disadvantage of the representation (|2.4j) is that it is valid only for 5fts > 1 . 
This can be improved by using a complex contour integral representation. Starting 
point is the following representation for the Gamma- function |46) . 



Lemma 2.6. For z 4 Z we have 



T(z) 



2i siniirz) 



(-t)* -1 e-*dt, 



where the anticlockwise contour C consists of a circle C3 of radius e < 2-k and 
straight lines C\, respectively Ci, just above, respectively just below, the x-axis; see 
Figure [IJ 
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Figure 1 . Contour C in Lemma 123 



Proof. Assume 5ftz > 1. As the integrand remains bounded along C3, no contribu- 
tions will result as e — > 0. Along C\ and C2 we parameterize as given in Figure Q] 
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and thus for 3?z > 1 

oo 

lim [(-ty-ie-'dt = f e.- t7T{z - 1) u z - l e- u du+ I e l < z -^ 



u z - x e- u 



du 



u z ~ 1 e~ u (e mz ~ e~ inz ) du 

oo 

= -2isin(7rz) / u z ~ l e~ u du, 
o 
which implies the assertion by analytical continuation. □ 

This representation for the Gamma-function can be used to show the following 
result for the Hurwitz zeta function. 

Theorem 2.7. For s € C, s £ N, we have 

_ . . r(l-a) f (-t)"- 1 e- ta , 

c 
with the contour C given in Figure [7J 

Proof. We follow the previous calculation to note 

oo 

/l_£\s-l -ta p -ta 
- — —dt = -2ism(irs) / /J 8 " 1 -dt, 
1 - e l J 1 - e J 
c o 

and we use [15] 

sin(7rs)r(s) = _ 

to conclude the assertion. □ 

From here, properties previously given can be easily derived. For s <E Z the 
integrand does not have a branch cut and the integral can easily be evaluated using 
the residue theorem. The only possible singularity enclosed is at t = and to read 
off the residue we use the expansion 

n=0 

2.2. Barnes zeta function. The Barnes zeta function is a multidimensional gen- 
eralization of the Hurwitz zeta function. 

Definition 2.8. Let s € C with 3?s > d and c e K + , r € R^. The Barnes zeta 
function is defined as 

(2-6) C B (*,c|r*)= J2 7 T^ ^ - 

' — ' c + 77i • rr 
meNj 

If c = it is understood that the summation ranges over tti =/= 0. 

For f = l c i '■= (1, 1, ..., 1, 1), the Barnes zeta function can be expanded in terms 
of the Hurwitz zeta function. 
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Example 2.9. Let us consider d = 2 and r = (1,1). Then 

-. ^ 1 _ v^ fc + 1 _ ^ fc + c + 1 — c 

CB(S ' C|2) " *V C + mi+m9) «<° W « (C + fc)S 
= Cfl(a-l,c) + (l-c)CH(a,c). 

Example 2.10. Let ej. be the number of possibilities to write an integer fc as a 
sum over d non-negative integers. We then can write 



oo 



Cb{s,c\U)= V T ~ ; ; rj = V ef , M » ' 

■^^ (c + mi + ... + tow r ^^ (c+fc) s 
The coefficient e k can be determined for example as follows. Consider 



1 1 1 / oo \ / oo 



(i-x) d i-x i-x [2^ x j'"[2-. 



(«9„.fc 



On the other side, using the binomial expansion 



1 ^ T(d + k) k _^ (d + k-l)\ k 

2^ r(H\k\ x ~1^ ' 



(l-x) d ^ T(d)k\ ^ (d-l)!fc! 

^ fd + k-l\ k 

= M d -i ) x - 

This shows 

'd + fc-l\ 1 



( B (s,c\l d ) = J2 



d-1 J(c + k)» 

which, once the dimension d is specified, allows to write the Barnes zeta function 
as a sum of Hurwitz zeta functions along the lines in Example [ 



It is possible to obtain similar formulas for n rational numbers |27[ 128] . 
For some properties of the Barnes zeta function the use of complex contour 
integral representations turns out to be the best strategy. 

Theorem 2.11. We have the following representations: 

oo 

1 f e~ ct 
C B (s,c\r) = — r / t 8 - 1 —; dt 



T(l-s) 



^-ct 



(-ty- l 7^r— —dt, 



2 ™ f ' n-=i(i-e- 

with the contour C given in Figure [7J 

Exercise 2. Use equation (|2.5p . respectively Lemma \2. 61 to proof Theorem [2. Ill 

The residues of the Barnes zeta function and its values at non-positive integers 
are best described using generalized Bernoulli polynomials [70] . 
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Definition 2.12. We define the generalized Bernoulli polynomials S„ (x\r) by the 
equation 

e -x* (-l) d ^ (_t)»-d 



m=i(i-e-^) n.i^^o 



^^ 71.! 



Using Definition 12.121 in Theorem 12.111 one immediately obtains the following 
properties of the Barnes zeta function. 

Theorem 2.13. We have 

w 

Res C B (z,c\r) = J'^ . B^ z (c\f), z = 1,2, ...,d, 



(2) 



«^~ iJ^r/ & m 



Exercise 3. Use the first representation of £s(s, c\r) in Thcorcm l2 . 1 1 1 together with 
Definition 12.121 to show Theorem 12.131 Follow the steps of the proof in Theorem 
1231 

Exercise 4. Use the second representation of Cs{s,c\r) in Theorem 12.111 together 
with Definition 12.121 and the residue theorem to show Theorem 12.131 

2.3. Epstein zeta function. We now consider zeta functions associated with sums 
of squares of integers [38l [39] . 

Definition 2.14. Let s e C with 5Rs > d/2 and ceR+,fe M+. The Epstein zeta 
function is defined as 

Ce(s,c\r)= } ^ 5 ;r— . 

fr^ 1 d (c + rimf + r 2 m% + ... + r d m^) s 

If c = it is understood that the summation ranges over to =/= 0. 
Lemma 2.15. For !fts > d/2, we have 

oo 



Proof. This follows as before from property (|2.5[) of the Gamma-function. D 

As we have noted in the proof of Theorem 12.51 it is the small-t behavior of the 
integrand that determines residues of the zeta function and special function values. 
The way the integrand is written in Lemma 12.151 this t — > behavior is not easily 
read off. A suitable representation is obtained by using the Poisson resummation 

Pi- 
Lemma 2.16. Let reC with 5Rr > and t G M+, then 

oo i oo 

£ *** - S E ■-*"■ 

I — — OO I — — OO 
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Exercise 5. If F(x) is continuous such that 

\F(x)\dx < oo, 
then we define its Fourier transform by 

oo 

Hu) - / n*- * 

— OO 

If 

/ |F(u)| dw < oo, 

— oo 

then we have the Fourier inversion formula 

oo 

""» - / *"> e """ du - 

— oo 

Show the following Theorem: Let F € L 1 (M). Suppose that the series 

converges absolutely and uniformly in v, and that 

E l^(m)|<oo. 

Then 

J2 F{n + v) = J2 F{n)e 2mv . 

Hint: Note that 

G(v) = VJ F(n + v) 

nGZ 

is a function of v of period 1. 

Exercise 6. Apply Exercise [5] with a suitable function F(x) to show the Poisson 
rcsummation formula Lcmma l2.16l 

In Lemma 12.161 it is clearly seen that the only term on the right hand side that 
is not exponentially damped as t — > comes from the I = term. Using the 
rcsummation formula for all d sums in Lemma 12.151 after resumming the rh = 
term contributes 

oo 



oo 
df2 r 

t s-d/2-l e -ct dt 



TT° 


'/2 




I"- 




Vn ■ ■ ■ 


'',/ 


r( s ) 


n d/2 




T{, 


"I) 



y/n ■ ■ ■ r d T(s)c s - d / 2 
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All other contributions after resummation are exponentially damped as t — > and 
can be given in terms of modified Besscl functions |46| . 

Definition 2.17. Let 5ftz 2 > 0, then we define the modified Bessel function K v {z) 
as 

oo 

o 

Performing the resummation in Lemma 12.151 according to Lemma 12. 161 with 
Definition ^. 17l one obtains the following representation of the Epstein zeta function 
valid in the whole complex plane [3J] [75] . 

Theorem 2.18. We have 

- d/2 r( S -f) # _ Sj 2n°c 



d-2s 



Q(s,c\r) = 



V r i • • • r <i r(s) r(s) 1 /n • • • r d 



n£Z d /{0} 



— + ... + — 
n r d 



K d 8 [2wV~c(^ + ... + ^ 
2 \ \n r d 



Exercise 7. Show Theorem 12. 181 along the lines indicated. 

From Definition 12. 171 it is clear that the Bessel function is exponentially damped 
for large 5ftz 2 . As a result the above representation is numerically very effective as 
long as the argument of K d / 2 - s is large. The terms involving the Bessel functions 
are analytic for all values of s, the first term contains poles. As an immediate 
consequence of the properties of the Gamma-function one can show the following 
properties of the Epstein zeta function. 

Theorem 2.19. For d even, (^£(s,c\r) has poles at s = 4, | — 1, ..., 1, whereas for 
d odd they are located ais=|,| — 1, ..., 4, — t > '£%. Furthermore, 




for d odd 
for d even. 

Exercise 8. Use Theorem 12.181 and properties of the Gamma- function to show 
Theorem [2H 

This concludes the list of examples for zeta functions to be considered in what 
follows. A natural question is what the motivations arc to consider these zeta 
functions. Before we describe a few aspects relating to this question let us mention 
how all these zeta functions, and many others, result from a common principle. 

3. Boundary value problems and associated zeta functions 

In this section wc explain how the considered zeta functions, and others, are all 
associated with eigenvalue problems of (partial) differential operators. 
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Example 3.1. Let M = [0, L] be some interval and consider the Dirichlet boundary 
value problem. 

d 2 

P<t>n{X) := —X-2<f>n(x) = \ n (f> n (x), 0„(O) = (j) n (L) = 0. 

The solutions to the boundary value problem have the general form 

4>n(x) = Asin(\/\ n x) + Bcos(^/X n x). 
Imposing the Dirichlet boundary condition shows we need 

n (O) = B = O, ^> n (L) = Asm(L^) = 0, 

which implies 

rr7r z 
A„ = -j2~, n € N. 

We only need to consider n € N because non-positive integers lead to linearly 
dependent eigenfunctions. The zcta function Cp( s ) associated with this boundary 
value problem is defined to be the sum over all eigenvalues raised to the power 
(— s), namely 



00 1 

Cp( S ) = ^A^ s , 5R S >-. 



n=l 

So here the associated zeta function is a multiple of the zeta function of Ricmann, 

00 _ / r \ 2s 

w) = E(t) = 7 <-(")■ 

n— 1 x ' 

Example 3.2. The previous example can be easily generalized to higher dimen- 
sions. We consider explicitly two dimensions; for the higher dimensional situation 
see pQ. Let M = {(x, y)\x <G [0,Li],y <G [0, L 2 ]}. We consider the boundary value 
problem with Dirichlet boundary conditions on M, that is 

f d 2 d 2 \ 
P(f>n,m{x,y) = I -q-2 - ^-j + c I 0„, m (x,y) = A„, m ^„, m (a;,2/), 

4>7i,m(0,y) = <f>n,m(Ll,y) = <fi n ,m (x,0) = 4>n,m(x,L 2 ) = 0. 

Using the process of separation of variables, eigenfunctions are seen to be 

{ nirx\ fmiry\ 
<Pn,m{x,y) =Asin — — I sin I — — I , 



with the eigenvalues 



= (l7) + ( W +c ' n ' meN - 



The associated zeta function therefore is 

oo oo 

( P (s) = E E 



n— 1 m— 1 



2 / \ 2 

LJ ' \L 2 ■ +C 
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which can be expressed in terms of the Epstein zeta function given in Definition 
12.141 as follows, 



Cp(s) 



1 



(3-1) 



-Xe s,c 



- 4 CeUc 



n \ I it 
L~i) '\M 



i 

r 



Example 3.3. Similarly one can consider periodic boundary conditions instead of 
Dirichlct boundary conditions, this means the manifold M is given by M = S 1 x S . 
In this case the eigcnfunctions have to satisfy 

d d 

0n,m(O,y) = (/) nim (Ll,y), -K-4>n, m (0,y) = Q-<l>n,m(Li,y), 

d d 

4>n : m(x,0) = <p n , m (x,L 2 ), TT 4>n,m [X, 0) = — (|) n . m (l, L 2 ) ■ 

dy dy 

This shows 

(j>n, m {x,y) = Ae~ x e~^ v , 

which implies for the eigenvalues 

/ 2im \ / 2-kvh \ 

Xnm= \—) + \— 2 
The associated zeta function therefore is 

Cp{s) = ( £ (s,c\f), r = 

Clearly, in d dimensions one finds 

(p(s) = Ce (s,c\r), r = 



■ 



+ c, (n, m) € 1 



27tY 



Example 3.4. As a final example we consider the Schrodinger equation of atoms in 
a harmonic oscillator potential. In this case M = R 3 , and the eigenvalue equation 
reads 



a ^^ / 2 2 2^ 

--— A+ — Uoxx +uj 2 y + uj 3 z J 
2ni 2 



i(X>yj z ) / ^n 1 ,ri2,n3 ( Pn 1 ,n2,n3 (X, J/, Z). 



This differential equation is augmented by the condition that eigcnfunctions must 
be square integrable, 4> ni ,n 2 ,n 3 (x,y, z) G £ 2 (M. 3 ). As is well known, this gives the 
eigenvalues 



An 



fuvi [ n\ + - j + /tW2 ( n 2 + - j + frw 3 f n 3 + - 



(ni,n 2 ,n 3 ) e N 3 ,. 



This clearly leads to the Barnes zeta function 

(p(s) =Cs(s,c|f), 

where 

1 



|fe(wi + W2 + W3), f = fi(Wl,W2,U>3) . 



If M = R is chosen the Hurwitz zeta function results. 
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The above examples illustrate how the zeta functions considered in Section 2 
are all related in a natural way to eigenvalues of specific boundary value problems. 
In fact, zeta functions in a much more general context are studied in great detail. 
For our purposes the relevant setting is the setting of Laplace-type operators on a 
Riemannian manifold M, possibly with a boundary dM. Laplace-type means the 
operator P can be written as 

P = -g jk Vjvl - E, 

where gi k is the metric of M, V v is the connection on M acting on a smooth 
vector bundle V over M, and where E is an endomorphism of V . Imposing suitable 
boundary conditions, eigenvalues A„ and eigenfunctions <p n do exist, 

P<j>n(x) = X n (j)n{x), 

and assuming A„ > the zeta function is defined to be 

oo 

( p (s) = J2k s 

n=l 

for 5Rs sufficiently large. If there are modes with A„ = those have to be excluded 
from the sum. Also, if finitely many eigenvalues are negative the zeta function 
can be defined by choosing nonstandard definitions of the principal value for the 
argument of complex numbers, but we will not need to consider those cases. 

4. (Some) Motivations to consider zeta functions 

There are many situations where properties of zeta functions in the above context 
of Laplace-type operators are needed. In the following we present a few of them, but 
many more can be found for example in the context of number theory [3[ [H [23l [79] 
and quantum field theory [El C31 [El HS1 HE1 IH01 ISU E3 SU S21 IS3 IZ3J - 

4.1. Can one hear the shape of a drum? Let M be a two-dimensional mem- 
brane representing a drum with boundary dM . The drum is fixed along its bound- 
ary. Then possible vibrations of the drum and its fundamental tones are described 
by the eigenvalue problem 

d 2 d 2 \ 

"dx 2 ~ ~dtf) ^ x,y > = A «w>y)> Mx,y)\ {Xiy)edM = 0. 

Here, (x,y) denotes the variables in the plane, the eigenfunctions 4> n (x,y) describe 
the amplitude of the vibrations and A„ its fundamental tones. In 1966 Kac |56j 
asked if just by listening with a perfect ear, so by knowing all the fundamental 
tones A„ , it is possible to hear the shape of the drum. One problem in answering 
this question is, of course, that in general it will be impossible to write down the 
eigenvalues A„ in a closed form and to read off relations with the shape of the drum 
directly. Instead one has to organize the spectrum intelligently in form of a spectral 
function to reveal relationships between the eigenvalues and the shape of the drum. 
In this context a particularly fruitful spectral function is the heat kernel 

n=l 

which as t tends to zero clearly diverges. Given that some relations between the 
fundamental tones and properties of the drum arc hidden in the t — >• behavior 
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let us consider this asymptotic behavior very closely. Before we come back to the 
setting of the drum, let us use a few examples to get an idea what the structure of 
the t — >• behavior of the heat kernel is expected to be. 

Example 4.1. Let M = S 1 be the circle with circumference L and let P = 
—d 2 /dx 2 . Imposing periodic boundary conditions eigenvalues are 

27rfc x 2 



\k= I —J , fcez, 
and the heat kernel reads 

oo 

K s ,f) = £ e-(^T\ 

k— — oo 

From Lemma 12.161 we find the t — >• behavior 

Kgi(t) = L + (exponentially damped terms) . 

v47ri 

Note that with obvious notation this could be written as 

Kgi(t) = vol ( M) + (exponentially damped terms) . 

V47rt 

Example 4.2. The heat kernel for the d-dimensional manifold M — S 1 x • • • x S 1 
with P = — A clearly gives a product of the above and thus 

Km if) — Kgi (t) x • • • x K s i (t) = -T-^-vol(M) + (exponentially damped terms) . 

Example 4.3. To avoid the impression that there is always just one term that is 
not exponentially damped consider M as above but P = — A + to 2 . Then 

K(t) = e~ m t KM{t) = e" m "* ( -^vol(M) + exponentially damped terms 

1 ~(_i) 



\{M) V^ — - — m 2e t e 2 + (exponentially damped terms) 



{4n)d /2 V ^ £] 

In fact, the structure of the heat kernel observed in this last example is the 
structure observed for the general class of Laplace-type operators. 

Theorem 4.4. Let M be a d- dimensional smooth compact Riemannian manifold 
without boundary and let 

P = -g ]k V]V v k - E, 

where gi is the metric of M, \7 V is the connection on M acting on a smooth vector 
bundle V over M, and where E is an endomorphism of V. Then as t — > 0, 



K(t)~Y. a k tk 



±k-d/2 

with the so-called heat kernel coefficients a^. 

Proof. See, e.g., [44]. □ 
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In Example 14.31 one sees that 



1 ("I)* 2k. 



ak = J^yn-kt m vol(M) ' 

In general, the heat kernel coefficients are significantly more complicated and they 
depend upon the geometry of the manifold M and the endomorphism E |44) . 

Up to this point we have only considered manifolds without boundary. In order to 
consider in more detail questions relating to the drum, let us now see what relevant 
changes in the structure of the small-i heat kernel expansion occur if boundaries 
are present. 

Example 4.5. Let M = [0, L] and P = —d 2 /dx 2 with Dirichlct boundary condi- 
tions imposed. Normalized eigenfunctions are then given by 

2 . firlx\ 



<Pi(x) = \j — sin 



7r€ X 2 



and the associated eigenvalues are 

Using Lemma \2. 161 this time we obtain 

(4.1) K(t) = vol(M) h (exponentially damped terms). 

V47rt 2 

Notice that in contrast to previous results we have integer and half-integer powers 
in t occurring. 

Exercise 9. There is a more general version of the Poisson resummation formula 
than the one given in Lemma 12.161 namely 



(4-2) E *- t(i+Cr - \J E e 



■ 

t 



C=-oo 



Apply Exercise [S] with a suitable function F(x) to show equation (J4.2I) . 

Exercise 10. Consider the setting described in Example 14. 5 1 The local heat kernel 
is defined as the solution of the equation 

d d 2 



with the initial condition 



at bj.k^v)- 



KmK(t,x,y) =S(x,y). 

In terms of the quantities introduced in Example 14.51 it can be written as 

oo 

K(t,x,y) = J2<pe(*)<Pt(y)e~ Xtt - 

Use the resummation (|4.2p for K(t, x, y) and the fact that 

L 


to rediscover the above result (|4.ip . 
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Exercise 11. Let M = [0,L] and 

d 2 2 

with Dirichlet boundary conditions imposed. Find the small-i asymptotics of the 
heat kernel. 

Exercise 12. Let M = [0, L] x S 1 x ■ • • x S 1 be a d-dimcnsional manifold and 

d 2 

Impose Dirichlet boundary conditions on [0,L] and periodic boundary conditions 
on the circle factors. Find the small-i asymptotics of the heat kernel. 

As the above examples and exercises suggest, one has the following result. 

Theorem 4.6. Let M be a d- dimensional smooth compact Riemannian manifold 
with smooth boundary and let 

P = -g ]k V]V v k - E, 

where g^ k is the metric of M, \7 V is the connection on M acting on a smooth vector 
bundle V over M, and where E is an endomorphism of V . We impose Dirichlet 
boundary conditions. Then as t — > 0, 

oo 

K{t) ~ J2 ak tk ~ d/2 

fe=0, 1/2,1,... 

with the heat kernel coefficients a^. 

Proof. See, e.g., 03 . □ 

As for the manifold without boundary case, Theorem 14.41 the heat kernel coef- 
ficients depend upon the geometry of the manifold M and the endomorphism E, 
and in addition on the geometry of the boundary. Note, however, that in contrast 
to Theorem 14.41 the small-t expansion contains integer and half- integer powers in t. 

The same structure of the small-t asymptotics is found for other boundary con- 
ditions like Neumann or Robin, see [44], and the coefficients then also depend on 
the boundary condition chosen. In particular, for Dirichlet boundary conditions 
one can show the identities 

(4.3) a n = (47r)- d / 2 vol(Af), a 1/2 - {Ak)-^' 2 (--\ vol(dM), 

a result going back to McKean and Singer [66] . In the context of the drum, what 
the formula shows is that by listening with a perfect ear one can indeed hear certain 
properties like the area of the drum and the circumference of its boundary. But as 
has been shown by Gordon, Webb and Wolpert [45], one cannot hear all details of 
the shape. 



Exercise 13. Use Exercise [12] to verify the general formulas (|4.3j) for the heat 
kernel coefficients. 
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Instead of using the heat kernel coefficients to make the above statements, one 
could equally well have used zeta function properties for equivalent statements. 
Consider the setting of Theorem 14.61 The associated zeta function is given by 



71 = 1 

where it follows from Weyl's law [80l [8T] that this series is convergent for Sfts > d/2. 
The zeta function is related with the heat kernel by 

oo 

(4.4) { P ( 8 ) = f ^rjl*- 1 K(t)dt, 

o 

where equation (|2.5[) has been used. This equation allows us to relate residues and 
function values at certain points with the small-i behavior of the heat kernel. In 
detail, 

tA*\ t> , t \ a i~ z d d-1 1 2n + l 

(4-5) Res C p(,) = ^y, z= -, — , ...,-,- — , n e N , 

(4.6) CpHz) = {-l) q q\a i+a , qeK 



■Jo- 



Keeping in mind the vanishing of the heat kernel coefficients a* with half-integer 
index for dM = 0, see Theorem 14.4) this means for d even the poles are actually 
located only at z = d/2, d/2 — 1, ..., 1. In addition, for d odd we get (p(—q) = for 

geN . 



Exercise 14. Use Theorem 14.61 and proceed along the lines indicated in the proof 
of Theorem 12.51 to show equations (|4.5j) and (|4.6j) . 



Going back to the setting of the drum properties of the zeta function relate with 
the geometry of the surface. In particular, from (|4.3|) and (|4.5j) one can show the 
identities 

. , N vol(M) /l\ vol(dM) 

Res Cp(1) = -^, Res Cp [^ = — ^, 

and the remarks below equation (|4.3p could be repeated. 

4.2. What does the Casimir effect know about a boundary? We next con- 
sider an application in the context of quantum field theory in finite systems. The 
importance of this topic lies in the fact that in recent years, progress in many fields 
has been triggered by the continuing miniaturization of all kinds of technical de- 
vices. As the separation between components of various systems tends towards the 
nanometer range, there is a growing need to understand every possible detail of 
quantum effects due to the small sizes involved. 

Very generally speaking, effects resulting from the finite extension of systems 
and from their precise form are known as the Casimir effect. In modern technical 
devices this effect is responsible for up to 10% of the forces encountered in micro- 
electromechanical systems [THl [SO] . Casimir forces are of direct practical relevance 
in nanotechnology where, e.g., sticking of mobile components in micromachincs 
might be caused by them [76]. Instead of fighting the occurrence of the effect in 
technological devices, the tendency is now to try and take technological advantage 
of the effect. 
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Experimental progress in recent years has been impressive and for some config- 
urations allows for a detailed comparison with theoretical predictions. The best 
tested situations are those of parallel plates [H] and of a plate and a sphere 
[20l [2T1 l62l [63l [69] ; recently also a plate and a cylinder has been considered [13j [37] . 
Experimental data and theoretical predictions are in excellent agreement, see, e.g., 
[8] I25 [ IB"T ] 164 ] . This interplay between theory and experiments, and the intriguing 
technological applications possible, are the main reasons for the heightened interest 
in this effect in recent years. 

In its original form, the effect refers to the situation of two uncharged, parallel, 
perfectly conducting plates. As predicted by Casimir [17], the plates should attract 
with a force per unit area, F(a) ~ 1/& 4 , where a is the distance between the plates. 
Two decades later Boyer [10] found a repulsive pressure of magnitude F(R) ~ 1/R 4 
for a perfectly conducting spherical shell of radius R. Up to this day an intuitive 
understanding of the opposite signs found is lacking. One of the main questions in 
the context of the Casimir effect therefore is how the occurring forces depend on the 
geometrical properties of the system considered. Said differently, the question is 
" What does the Casimir effect know about a boundary?" In the absence of general 
answers one approach consists in accumulating further knowledge by adding bits 
of understanding based on specific calculations for specific configurations. Several 
examples will be provided in this section and we will see the dominant role the 
zeta functions introduced play. However, before we come to specific settings let us 
briefly introduce the zeta function regularization of the Casimir energy and force 
that we will use later on. 

We will consider the Casimir effect in a quantum field theory of a non-interacting 
scalar field under external conditions. The action in this case is |55j 

(4.7) S[&] =~l f *(&) (A - V(x)) $(x) dx 

M 

describing a scalar field $(x) in the background potential V(x). We assume the 
Riemannian manifold M to be of the form M = S 1 x M s , where the circle S 1 
of radius (3 is used to describe finite temperature T — 1/(3 and M s , in general, 
is a d-dimensional Riemannian manifold with boundary. For the action (|4.7[) the 
corresponding field equations are 

(4.8) (A - V(x))$(x) = 0. 

If M s has a boundary dM s , these equations of motion have to be supplemented by 
boundary conditions on dM s . Along the circle, for a scalar field, periodic boundary 
conditions are imposed. 

Physical properties like the Casimir energy of the system are conveniently de- 
scribed by means of the path-integral functionals 

(4.9) Z[V]= f e- sm D$, 

where we have neglected an infinite normalization constant, and the functional 
integral is to be taken over all fields satisfying the boundary conditions. Formally, 
equation (|4.9[) is easily evaluated to be 



(4.10) T[V] = -\nZ[V] = ^\ndct[(-A + V(x))/fi 2 
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where [i is an arbitrary parameter with dimension of a mass to adjust the dimension 
of the arguments of the logarithm. 

Exercise 15. In order to motivate equation (|4.10[) show that for P a positive 
definite Hcrmitian (N x iV)-matrix one has 

-k(*,P*)(dx) ^(dctP)- 1 / 2 , 

where (dx) = d n x(2ii)~ n l 2 . For P = —A + V(x) and interpreting the scalar product 
(x,Px) as an L 2 (M)-product, one is led to (|4.10p by identifying Z?$ with (dx). 

Clearly equation (|4.10[) is purely formal because the eigenvalues A„ of — A + 
V(x) grow without bound for n — > oo and thus expression (|4.10p needs further 
explanations. 

In order to motivate the basic definition let P be a Hcrmitian (N x iV)-matrix 
with positive eigenvalues A„. Clearly 



N , N 

d 



lndctP = ^]lnA rl = --r^2K S 






s=0 



and the determinant of P can be expressed in terms of the zeta function associated 
with P. This very same definition, namely 

(4.11) lndctP= -Cp(0) 

with 



oo 



(4-12) Cp(-s) = ^ A" s 

n=l 

is now applied to differential operators as in (|4.10[) . Here, the series representation 
is valid for 5fts large enough, and in (]4.11[) the unique analytical continuation of the 
series to a neighborhood about s = is used. 

This definition was first used by the mathematicians Ray and Singer |73j to 
give a definition of the Reidemeister-Franz torsion. In physics, this regularization 
scheme took its origin in ambiguities of dimensional regularization when applied 
to quantum field theory in curved spacetime [291 15 1J . For applications beyond the 
ones presented here see, e.g., [HI [15l [26l [30l EH HH HH [74] . 

The quantity T[V] is called the effective action and the argument V indicates 
the dependence of the effective action on the external fields. The Casimir energy is 
obtained from the effective action via 

9/3 L J 2 9/3 

Here, we will only consider the zero temperature Casimir energy 

(4.14) E Cas = hm E 

p— foo 

and we will next derive a suitable representation for Ecas ■ We want to concentrate 
on the influence of boundary conditions and therefore we set V(x) = 0. The relevant 
operator to be considered therefore is 

d 2 



( 4 -!3) E = —,nv}= — ^Cr,„>(0). 
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where r £ S 1 is the imaginary time and A s is the Laplace operator on M s . In 
order to analyze the zeta function associated with P we note that eigenfunctions, 
respectively eigenvalues, are of the form 

1 2^ 

<t>n,j(T,y) = -je 3 T ifj{y), 

with 

-KVj{y) =E?ipj(y), 
where y S M s . For the non-selfinteracting case considered here, Ej are the one- 
particle energy eigenvalues of the system. The relevant zeta function therefore has 
the structure 

(4.15) fr(.)_ tf.((?f) 2 + E>Y'. 

n=-oo j=l \ V H ' / 

We repeat the analysis outlined previously, namely we use equation (|2.5[) and we 
apply Lemma 12.161 to the n-summation. In this process the zeta function 



and the heat kernel 



oo 

-2s 
3 



J=l fc=0,l/2,l,... 

of the spatial section are the most natural quantities to represent the answer, 

oo 

^ = W7) E t^e-i^tRpMdt 

/3 r( s -i/2) 

Cp s (s- 1/2) 



4tt F(s) 
J " 

7TF( S ) 

v ' n— 1 



„2.j2 



^ t s -V 2 e-^K Ps {t)dt. 



o 
For the Casimir energy we need (I? = d + 1) 

Cp/^0) = Cp(0) + Cp(0)ln M 2 

= -j3 {FP (p s (-1/2) + 2(1 - ln2)Rcs ( Ps (-1/2) 

OO 

n=l 

= -0 (Vp Cp.(-V2) - -4=ttfl/a [(In m 2 ) + 2(1 -In 2); 



(4.16) 



47T 



+ ^Ly lt-*' 2 e-{^)K Ps {t)dt, 



20 KLAUS KIRSTEN 

with the finite part FP of the zeta function and where equations (|4.5[) and (|4.6|) 
together with the fact that 

K M {t)=K s i(t) K Ps {t) 

have been used, in particular 

(4.17) Rc s<p ,(-I) = -^. fr<P)- -£-„/,. 

At T = we obtain for the Casimir energy, see equations (|4.13[) and (|4.14[) . 

(4.18) E Cas = hm E = ]-FP Cp.(-1/2) - -±=a D/2 lnjl 2 , 

with the scale fi = (fie/ 2). Equation (|4.18[) implies that as long as ao/2 t^ the 
Casimir energy contains a finite ambiguity and renormalization issues need to be 
discussed. Note from (|4.17|) that whenever (p s (— 1/2) is finite no ambiguity exists 
because cld/2 — 0. In the specific examples chosen later we will make sure that 
these ambiguities are absent and therefore a discussion of renormalization will be 
unnecessary. 

In a purely formal calculation one essentially is also led to equation (|4.18|) . As 
mentioned, in the quantum field theory of a free scalar field the eigenvalues of a 
Laplacian are the square of the energies of the quantum fluctuations. Writing the 
Casimir energy as (one-half) the sum over the energy of all quantum fluctuations 
one has 



-. oo 

(4-19) Eca^-Y,* 



1/2 

2^'" k 



k=0 



and a formal identification 'shows' 



(4-20) Ecas = \tP. (-\ 

Clearly, the expression (|4.19[) is purely formal as the series diverges. However, when 
(p B (— 1/2) turns out to be finite this formal identification yields the correct result. 
Otherwise, the ambiguities given in (|4.18p remain as discussed above. 

An alternative discussion leading to definition (|4.18[) can be found in [7]. 

As a first example let us consider the configuration of two parallel plates a dis- 
tance a apart analyzed originally by Casimir |17j . For simplicity we concentrate on 
a scalar field instead of the electromagnetic field and we impose Dirichlet boundary 
conditions on the plates. The boundary value problem to be solved therefore is 

-Av,k(x, y, z) = \ku k (x, y, z) 

with 

u>k(0,y,z) = u k (a,y,z) = 0. 

For the time being, we compactify the (y, z)-directions to a torus with perimeter 
length R and impose periodic boundary conditions in these directions. Later on, 
the limit R — > oo is performed to recover the parallel plate configuration. Using 
separation of variables one obtains normalized eigenfunctions in the form 



u ei i 2 i(x,y,z) = \l~/^2 sin [— x 
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with eigenvalues 



A/ 



/2tt£i 



2-kL 



V R J \ R 
This means we have to study the zeta function 

(4.2i) cw= E E 



— , (4,^ 2 ) eZ", ^e 



«=1 L 



2tt^i 



i? 



2tt^ 



As i? — > 00 the Ricmann sum turns into an integral and we compute using polar 
coordinates in the (y, z ) -plane 



CM 



9tt / 2^ 



00 00 r 



— 00 — 00 

00 



I- 



dfe dfci 



s)S> 



1=1 



k 2 + - 



dk 



R 2 1 




2tt 2(1- 
R 2 


*) 


4tt(1- 
i? 2 


*) 



E 
E 



nf 
a 

2(-s+l) 



k 2 + — 



a 

2-2s 



-s+1 



C«(2s - 2). 



47r(l — sj \a/ 
T/2 as needed for the Casimir energy we obtain 
1\ & 2 /tt\3 , , „ s JJV 

2y ~ 4?r 3 \o- 
The resulting Casimir force per area is 



Setting s 
(4.22) 



( 



er<«<-3) 



720a 3 



(4.23) 



-Fbas 



9 £ C a S 



9a i? 2 480a 4 ' 

Note, that this computation takes into account only those quantum fluctuations 
from between the plates. But in order to find the force acting on the, say, right 
plate the contribution from the right to this plate also has to be counted. To find 
this part we place another plate at the position x — L where at the end we take 
L — > 00. Following the above calculation, we simply have to replace a by L — a to 
see that the associated zeta function produces 

^ \ 2) 720(L-a) 3 

and the contribution to the force on the plate at x = a reads 

7T 2 

Cas ~ 480(L-a) 4 ' 

This shows the plate at x = a is always attracted to the closer plate. As L — > 00 it 
is seen that equation (|4.23[) also describes the total force on the plate at x = a for 
the parallel plate configuration. 
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Exercise 16. Consider the Casimir energy that results in the previous discussion 
when the compactification length R is kept finite. Use Lemma \2. 181 to give closed 
answers for the energy and the resulting force. Can the force change sign depending 
on a and R? 

More realistically plates will have a finite extension. An interesting setting that 
we are able to analyze with the tools provided are pistons. These have received 
an increasing amount of interest because they allow the unambiguous prediction of 

forces nn na UHl EH E3 ■ 

Instead of having parallel plates let us consider a box with side lengths L\,L 2 
and L 3 . Although it is possible to find the Casimir force acting on the plate at 
x = L\ resulting from the interior of the box, the exterior problem has remained 
unsolved until today. No analytical procedure is known that allows to obtain the 
Casimir energy or force for the outside of the box. This problem is avoided by 
adding on another box with side lengths L — L\,L 2 and L 3 such that the wall at 
x = L\ subdivides the bigger box into two chambers. The wall at x = L\ is assumed 
to be movable and is called the piston. Each chamber can be dealt with separately 
and total energies and forces are obtained by adding up the two contributions. 
Assuming again Dirichlct boundary conditions and starting with the left chamber, 
the relevant spectrum reads 



(4.24) X e 



■nix 
£1 



7r£2 

L 2 



and the associated zeta function is 



(4.25) 



C(«) 



E 



«i,4AM L 



£1 



7^3 
L 3 

L 2 



h,l 2 ,h GN, 



7T4 

L 3 



One way to proceed is to rewrite (|4.25[) in terms of the Epstein zeta function in 
Definition [HIl 



Exercise 17. Use Lemma [2.181 in order to find the Casimir energy for the inside 
of the box with side lengths L\, L 2 and L 3 and with Dirichlct boundary conditions 
imposed. 

Instead of using Lemma [2.181 we proceed as follows. We write first 



COO 



5 E E 



(4.26) 



l=-00 <!,<3=1 
TTi2 

L 2 






1 OO 



7T^3 

L 3 



This shows that it is convenient to introduce 

2 



(4.27) 



&w= E 



1=1 



7T4 

L 2 



L 2 



^3 
L 3 



7T^3 
£3 



We note that this could be expressed in terms of the Epstein zeta function given 
in Definition 12.141 However, it will turn out that this is unnecessary. 
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Also, to simplify the notation let us introduce 
Using equation (|2.5|) for the first line in (|4.26p we continue 



2 
Mats 



v ; « 1= -oo t 2 ,t 3 =± 



t s x exp ■ 



-f 



71^ 



Vizis 



(It 



We now apply the Poisson resummation in Lemma 12.161 to the ^i-summation and 
therefore we get 



COO 



II 



20F r(s) 



E E 



t s 2 exp 



2/2 



L{i{ 



~ *#Lfe, *>■ dt 



(4.28) 



'&(*)■ 



The £i = term gives a Cc-t erm > the ^i ^0 terms are rewritten using (|2.17|) . The 
outcome reads 



(4.29) COO 



^-^ / i 

20F T(s) w V 2 



MH 



+ 



2£T! v 1 £ ' 

/5F r(s) H \iA , 

y ' £l,£ 2 ,£ 3 = l V *2*3 






Ky_ s {2L 1 t 1 m 2 i 3 ). 



We need the zeta function about s = —1/2 in order to find the Casimir energy and 
Casimir force. 

Let s = —1/2 + e. In order to expand equation (|4.29[) about e = we need to 
know the pole structure of Cc(s). From equation (|2.18|) it is expected that Cc( s ) 
has at most a first order pole at s = — 1/2 and that it is analytic about s — — 1. So 
for now let us simply assume the structure 



4 + £ 



= - e ^Cc(-l) + FPCc(-\ 



O(e), 



Ce(-l + e) = (c(-l) + e(i(-l)+0(e 2 ), 



where Res (c( — 1/2) and FP Cc(~ 1/2) will be determined later. With this structure 
assumed, we find 



( 



1 



(4.30) 



+§ (4I- 1 ) + &(-l)0"4 - 1)) - jFP & 



^ OO 



A^2*3 



K x {2L x t xiH ^) . 



il,£ 2 ,i3 = l 

This shows that the Casimir energy for this setting is unambiguously defined only 
if Ce(-l) = and Res Ce(-l/2) = 0. 
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Exercise 18. Show the following analytical continuation for Cc(s): 



+- 



2Lt 



-1/2 



rif r(s) 



EE 



^3 
7T& 



8-1/2 



A"i 



/2ttL 2 



V ^3 



2 = 1^3=1 

Read off that Cc(-l) = Res Cc(-l/2) = 0. 

Using the results from Exercise IT51 the Casimir energy, from equation (|4.30|l . can 
be expressed as 



(4.32) 



Ecu 



8tt cV ; 4 SC V 2 



-. OO 

2^ ?-< 



fl,4a,*3=l 



Wjfc 



«1 



K^LA^is). 



Exercise 19. Use representation (|4.3ip to give an explicit representation of the 
Casimir energy (|4.32|) . 



For the force this shows 



(4.33) F< 



Cas 



i i °° 



fl ^3^3=1 



Mf2<?3 



9Li 



KtfL&HM 



Exercise 20. Use Dchnition l2.17l to show that K u {x) is a monotonically decreasing 
function for x <E R+. 

Exercise 21. Determine the sign of Cc( — !)■ What is the sign of the Casimir force 
as L\ — > oo? What about L\ — > 0? 

Remember that the results given describe the contributions from the interior of 
the box only. The contributions from the right chamber are obtained by replacing 
L\ with L — L\. This shows for the right chamber 



Eca 



L-L 



^(-U-iFP&(-i 



-. CXD 

— y 

2tt ^ 



^e 2 £3 



£i 



K 1 (2(L-L l )hm 2 i 3 ), 



F Ca 



1 1 °° 



Wa4> 



*1 



9Li 



K"i(2(Z,-Li)^ a / 3 ). 



Adding up, the total force on the piston is 






2- 



(4.34) 



- y 

_. oo 

— y 



lAA=l 
OO 



Mf 2 f 3 



9Li 



#i(2Li*i Wa < 3 ) 



»,=i 



fj.e 2 e 3 



h 



—Kx^L-Lx^m^). 



This shows, using the results of Exercise [20l that the piston is always attracted to 
the closer wall. 
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Although we have presented the analysis for a piston with rectangular cross- 
section, our result in fact holds in much greater generality. The fact that we an- 
alyzed a rectangular cross-section manifests itself in the spectrum (|4.24[) . namely 
the part 

L 2 J \L 3 

is a direct consequence of it. If instead we had considered an arbitrary cross-section 
C, the relevant spectrum had the form 



where, assuming still Dirichlet boundary conditions on the boundary of the cross- 
section C, /i 2 is determined from 



d 2 d 2 



= 0. 
(v,z)edC 



Proceeding in the same way as before, replacing ^£ 2 e 3 with ^ and introducing Cc(s) 
as the zeta function for the cross-section, 

oo 

Cc(s)=X>- 2 *, 

equation (|4.28|l remains valid, as well as equations (|4.29[) and (|4.30|) . So also for an 
arbitrary cross-section the total force on the piston is described by equation (|4.34[) 
with the replacements given and the piston is attracted to the closest wall. 

Exercise 22. In going from equation (J4.28I) to (|4.29| l the fact that /i 2 e > is 
used. Above we used fj, 2 > which is true because we imposed Dirichlet boundary 
conditions. Modify the calculation if boundary conditions are chosen (like Neumann 
boundary conditions) that allow for do zero modes fj, 2 = [S5] ■ 

We have presented the piston set-up for three spatial dimensions, but a similar 
analysis can be performed in the presence of extra dimensions |58j . Once this 
kind of calculation is fully understood for the electromagnetic field it is hoped that 
future high-precision measurements of Casimir forces for simple configurations such 
as parallel plates can serve as a window into properties of the dimensions of the 
universe that are somewhat hidden from direct observations. 

As we have seen for the example of the piston, there are cases where an un- 
ambiguous prediction of Casimir forces is possible. Of course the set-up we have 
chosen was relatively simple and for many other configurations even the sign of 
Casimir forces is unknown. This is a very active field of research; some references 
are [5J [3H H31 ISZ1 GUI US] . Further discussion is provided in the Conclusions. 

5. BOSE-ElNSTEIN CONDENSATION OF BOSE GASES IN TRAPS 

We now turn to applications in statistical mechanics. We have chosen to ap- 
ply the techniques in a quantum mechanical system described by the Schrodingcr 
equation 

(5-1) f -— A + V(x,y,z) J 4> k (x,y,z) = X k 4>k(x,y,z), 
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that is we consider a gas of quantum particles of mass m under the influence of 
the potential V(x,y,z). Specifically, later we will consider in detail the harmonic 
oscillator potential 

V(x, y,z) = — (ujix 2 + uj 2 y 2 + w 3 z 2 ) 



briefly mentioned in Example 13.41 as well as a gas confined in a finite cavity. 

Thcrmodynamical properties of a bose gas, which is what we shall consider in 
the following, arc described by the (grand canonical) partition sum 

00 
(5.2) g = -^ln(l-e-^ (Afc -' j) 

A:=0 

where (3 is the inverse temperature and fi is the chemical potential. We assume 
the index k = labels the unique ground state, that is, the state with smallest 
energy eigenvalue Ao- From this partition sum all thermodynamical properties are 
obtained. For example the particle number is 



<«> "41 



OO 1 



^e/J(A fc -M)_i' 

T ' v k=0 



where the notation (dq/dfj,\T,v) indicates that the derivative has to be taken with 
temperature T and volume V kept fixed. The particle number is the most impor- 
tant quantity for the phenomenon of Bose-Einstein condensation. Although this 
phenomenon was predicted more than 80 years ago [9] [32] it was only relatively 
recently experimentally verified (2] [TTJ [24] . Bose-Einstein condensation is one of 
the most interesting properties of a system of bosons. Namely, under certain condi- 
tions it is possible to have a phase transition at a critical value of the temperature 
in which all of the bosons can condense into the ground state. In order to un- 
derstand at which temperature the phenomenon occurs a detailed study of N, or 
alternatively q, is warranted. This is the subject of this section. 

Wc first note that from the fact that the particle number in each state has to be 
non-negative it is clear that [i < Ao has to be imposed. It is seen in (|5.2[) that as 
(3 — > (high temperature limit) the behavior of q cannot be easily understood. But 
contour integral techniques together with the zeta function information provided 
makes the analysis feasible and it will allow for the determination of the critical 
temperature of the bose gas. 

Let us start by noting that from 



\n(l-x) = -J2—, iov\x\<l, 

n=l 

the partition sum can be rewritten as 

(5.4) g = ffl e -ft>. 

ra=l k=0 

The /3 — > behavior is best found using the following representation of the expo- 
nential. 

Exercise 23. Given that 

lim \T(x + iy)\ e^ |y| \yfi~ x = V~2^, x,yeE., 

|»|-»-oo 
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and 

T(z) = ^e( 2 -^) i°g- -* (i + (i)) i 

as \z\ — k oo, show that 

<x-Hoo 

(5.5) e~ a = — / a"* T(t) dt, 

2m J 

a—ioo 

valid for a > 0, |arg a| <§-<$, 0<<5< tt/2. 

Before we apply this result to the partition sum (|5.4[) let us use a simple example 
to show how this formula allows us to determine asymptotic behavior of certain 
series in a relatively straightforward fashion. From Lemma 12.161 we know that 



OO 

1=1 


1 oo 

= \ E ^ 

fc-oo 


1 

~ 2 




i rw i 

~ 2V^"2 + ^ 


h 



(5 - 6 > - iV?-5 + V^*" v - 

As j3 — > it is clear that the series on the left diverges and Lcmma ^. 161 shows that 
the leading behavior is described by a 1/v 7 ? term, followed by a constant term, 
followed by exponentially damped corrections. Let us see how we can easily find 
the polynomial behavior as (3 — > from (J5.5J1 . We first write 



(T + ZOO 



E^ £2 = Ei / w 2 rnt)dt. 

B 1 B 1 ^ J 



1=1 i=1 a- l0 o 



Here, a > is assumed by Exercise [23] However, in order to be allowed to inter- 
change summation and integration we need to impose a > 1/2 and find 



CT+ZOO 



E^ 2 - h r*r(*)5:<-** 

/3- t r(*)Cfl(2t)dt 



fcl CT -ioo fcl 



er-Hoo 

= J- /" a-' 

27rj 

c 

In order to find the small-/? behavior, the strategy now is to shift the contour to 
the left. In doing so we cross over poles of the integrand generating polynomial 
contributions in (3. For this example, the right most pole is at t = 1/2 (pole of the 
zeta function of Riemann) and the next pole is at t = (from the Gamma- function) . 
Those are all singularities present as (r(— 2n) = for n € N. Therefore, 

CT+ioo 



(=1 



CT+ZOO 
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where a < and where contributions from the horizontal lines between a ± ioo and 
a ± ioo are neglected. For the remaining contour integral plus the neglected hori- 
zontal lines one can actually show that they will produce the exponentially damped 
terms as given in (|5.6[) . How exactly this actually happens has been described in 
detail in |35j . 

Exercise 24. Argue how Y^=\ e _ ^™° , /3 > 0, a > 0, behaves as /? — > by using 
the above procedure. Determine the leading three terms in the expansion assuming 
that the contributions from the contour at infinity can be neglected. 

Exercise 25. Find the leading three terms of the small-/? behavior of 

£log(l-e-^) 

n=l 

assuming that the contributions from the contour at infinity can be neglected. 



Let us next apply the above ideas to the partition sum (|5.4[) . As a further 
warmup, for simplicity, let us first set \x = 0. Not specifying A^ for now and using 



fc=0 

for -fts > M large enough to make this series convergent, we write 






cr+ioc 




n=l k=Q 
OO OO -. -. ,. 

= EE-^" y (^n)-*r W * 

(T — ZOO 

cr+'i oo 

= in J ^"^WChC* + i)C(*)*. 

CT — 200 

Here u > M is needed for the algebraic manipulations to be allowed. It is clearly 
seen that the integrand has a double pole at t = 0. The right most pole (at M) 
therefore comes from £(£), and the location of this pole determines the leading 
/3 — > behavior of the partition sum. 

For the harmonic oscillator potential, in the notation of Example l3.4[ the Barnes 
zeta function occurs and we have 

(5.7) q= ^Ti J ^* r WCfl(i + l)CB(i,c|r)di. 

(7 — ioo 

The location of the poles and its residues are known for the Barnes zeta function, 
see Definition 12.121 and Theorem 12.131 in particular one has 

1 



Res Cg(3,c|r) = 



2frW ' 
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where, as is common, the geometric mean of the oscillator frequencies 

n = (wiw 2 w 3 ) 1/3 

has been used. The leading order of the partition sum therefore is 

Exercise 26. Use Definition 12.121 and Theorem 12.131 to find the sublcading order 
of the small-/3 expansion of the partition sum q. 

Exercise 27. Consider the harmonic oscillator potential in d dimensions and find 
the leading and sublcading order of the small-/? expansion of the partition sum q. 

If instead of considering a bosc gas in a trap we consider the gas in a finite three- 
dimensional cavity M with boundary dM we have to augment the Schrddingcr 
equation (|5.1I) by boundary conditions. We choose Dirichlet boundary conditions 
and thus the results for the heat kernel coefficients (|4.3|) arc valid. 

From equation (|4.5|) we conclude furthermore that the rightmost pole of £(s) is 
located at s = 3/2 and that 

RcsC<' 3 ^- a ° - VOl(M) 



2) r(i)- 4** ' 

furthermore the next pole is located at s = 1. For this case, the leading order of 
the partition sum therefore is 

? = ^C*(0voi(M) + O(/ ri). 

One way to read this result is that the bose gas does know the volume of its container 
because it can be found from the partition sum. This is completely analogous to 
the statement for the drum where we used the heat kernel instead of the partition 
sum. 

Subleading orders of the partition sum reveal more information about the cavity, 
see the following exercise. But as for the drums, the gas does not know all the details 
of the shape of the cavity because there are different cavities leading to the same 
eigenvalue spectrum [45] . Those cavities cannot be distinguished by the above 
analysis. 

Exercise 28. Consider a bose gas in a d-dimensional cavity M with boundary dM. 
Use (|4.3|) and (|4.5|) to find the leading and sublcading order of the small- j3 expansion 
of the partition sum q. What does the bose gas know about its container, meaning 
what information about the container can be read of from the high-temperature 
behavior of the partition sum? 

In order to examine the phenomenon of Bosc-Einstcin condensation we have 
to consider non-vanishing chemical potential. Close to the phase transition, as 
we will see, more and more particles have to reside in the ground state and the 
value of the chemical potential will be close to the smallest eigenvalue, which is the 
'critical' value for the chemical potential, /i c = Ao- Near the phase transition, for 
the expansion to be established, it will turn out advantageous to rewrite Afc — /i 
such that the small quantity fi c — fi appears, 

Afc - n = Afc - n c + He- (i = A fe - Ao + He - M- 
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Given the special role of the ground state, we separate off its contribution and write 



OO OO -j 

q = <Zo + T T - e-^^-^ e^"^""). 

n=l fc=l 



Note that the fc-sum starts with k = 1, which means that the ground state is not 
included in this summation. Employing the representation (|5.5|) only to the first 
exponential factor and proceeding as before we obtain 



Co(t)dt, 
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fc=l 

In order to determine the small- /3 behavior of expression (|5.8[) let us discuss the pole 
structure of the integrand. Given fi c — /i > 0, the polylogarithm Lii+^e - ^^ - ^) 
does not generate any poles. Concentrating on the harmonic oscillator, we find 

ResCoO) = ^3, 

1/1 1 1 
Res Co 2 = —2 + + 

2n z \LO1OJ2 W1W3 W2W3 

Note that Co( s ) i s the Barnes zeta function as given in Definition 12.81 with c = 
where we have to exclude rn = from the summation. However, clearly the residues 
at s = 3 and s = 2 can still be obtained from Theorem 12.131 with c — > taken. 
Shifting the contour to the left we now find 

+ ^ Lia ( e -^- M) ) (— + — + — ) + ■•• 

Z(ph) z \ / \W1CJ2 W1W3 OJ2UJ3 J 

In order to find the particle number N we need the relation for the polylogarithm 

dU n (x) 1 

which follows from (15.91). So 



" = N0+ (Bm Lh (^ { ^) 



1 rLi2 r e -«"-" > w 



2(/3h) 2 V / \w1cj2 ^1^3 W2W3 
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Exercise 29. Use ([575]) and (|5T9|1 to show 

Li„ (e _a: ) = 0?» - <«(« - 1) + ••■ 

valid for n > 2. What does the subleading term look like for n = 2? 

As the critical temperature is approached \l — > (x c and with Exercise [2H1 the 
particle number close to the transition temperature becomes 

< 5 - 10 ' N - ™" + 7i§ + l^(— + — + - 

(pfiiz)' 3 2{pn) z \cu1UJ2 W1W3 W2<JJ3 / 

The second and third term give the number of particles in the excited levels (at 
high temperature close to the phase transition). 

The critical temperature is defined as the temperature where all excited levels are 
completely filled such that lowering the temperature the ground state population 
will start to build up. This means the defining equation for the critical temperature 
T c = 1//9 C in the approximation considered is 



(/3 C M7) 3 ' 2(/3 c h) 2 \Ld1uj2 W1CJ3 W2W3 

Solving for /3 C one finds 

r-^-ssf*'*- 1 *}- 

Here, To is the critical temperature in the bulk limit (N — > 00) 

N \V3 

T = Ml 



,Cr(3) 
and 



2 \W1W2 0J1UI3 UJ2OJ2, , 

Different approaches can be used to obtain the same answers [47l ESI HH [50] . 

If only a few thousand particles are used in the experiment the finite- N correc- 
tion is actually quite important. For example the first successful experiments on 
Bose-Einstein condensates were done with rubidium [5] at frequencies ui\ — cj 2 = 
240tt/V8 s^ 1 and w 3 = 2407TS- 1 . With N = 2000 one finds T c - 31.9nK= 0.93 T 
[59] . a significant correction compared to the thcrmodynamical limit. 

Exercise 30. Consider the bose gas in a d-dimensional cavity. Find the particle 
number and the critical temperature along the lines described for the harmonic 
oscillator. What is the correction to the critical temperature caused by the finite 
size of the cavity? (For a solution to this problem see [50].) 

6. Conclusions 

In these lectures some basic zeta functions are introduced and used to analyze the 
Casimir effect and Bose-Einstein condensation for particular situations. The basic 
zeta functions considered are the Hurwitz, the Barnes and the Epstein zeta function. 
Although these zeta functions differ from each other they have one property in 
common: they are based upon a sequence of numbers that is explicitly known and 
given in closed form. The analysis of these zeta functions and of the indicated 
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applications in physics is heavily based on this explicit knowledge in that well- 
known summation formulas are used. 

In most cases, however, an explicit knowledge of the eigenvalues of, say, a Lapla- 
cian will not be available and an analysis of the associated zeta functions will be 
more complicated. In recent years a new class of examples where eigenvalues are 
defined implicitly as solutions to transcendental equations has become accessible. 
In some detail let us assume that eigenvalues are determined by equations of the 
form 

(6.1) Ftfa, n ) = 

with £, n suitable indices. For example when trying to find eigenvalues and eigen- 
functions of the Laplacian whenever possible one resorts to separation of variables 
and I and n would be suitable 'quantum numbers' labeling cigcnfunctions. To be 
specific consider a scalar field in a three dimensional ball of radius R with Dirichlct 
boundary conditions. The eigenvalues Afc for this situation, with A; as a multiindex, 
are thus determined through 

-A(f> k (x) = \ k <j)k{x), cj)(x)\\ x \ =B = 0. 

In terms of spherical coordinates (r, f2), a complete set of eigenfunctions may be 
given in the form 

where Yi, m (f2) arc spherical surface harmonics [40], and J v are Bessel functions 
of the first kind [46]. Eigenvalues of the Laplacian are determined as zeroes of 
Bessel functions. In particular, for a given angular momentum quantum number I, 
imposing Dirichlct boundary conditions, eigenvalues X^ n are determined by 

(6.2) J l+1/2 

Although some properties of the zeroes of Bessel functions are well understood [46] , 
there is no closed form for them available and we encounter the situation described 
by (|6.1j) . In order to find properties of the zeta function associated with this kind 
of boundary value problems the idea is to use the argument principle or Cauchy's 
residue theorem. For the situation of the ball one writes the zeta function in the 
form 



(6.3) C ( S ) = ^( 2 ; + l)_ / k- 2s — In, J l+1/2 (kR)dk, 

i=o m J y 

where the contour 7 runs counterclockwise and must enclose all solutions of (|6.2|) . 
The factor (2l + 1) represents the degeneracy for each angular momentum I and the 
summation is over all angular momenta. The integrand has singularities exactly at 
the eigenvalues and one can show that the residues are one such that the definition 
of the zeta function is recovered. More generally, in other coordinate systems, one 
would have, somewhat symbolically, 

(6.4) t(s) =E d ^ / k- 2s ^\nF 3 {k)dk, 

the task being to construct the analytical continuation of this object. The details 
of the procedure will depend very much on the properties of the special function Fj 
that enters, but often all the information needed can be found [57]. Nevertheless, 
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for many separable coordinate systems this program has not been performed but 
efforts are being made in order to obtain yet unknown precise values for the Casimir 
energy for various geometries. 
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